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Abstract. We present a framework for the dynamics of causal sets and coupled matter fields, 
which is a simplification and generalization of an approach we recently proposed. Given a set 
of fields including the gravitational one, the main step in implementing our proposal consists in 
writing their continuum-based action using as variables for the spacetime geometry the causal 
order and volume element. One then discretizes the resulting expression, with a procedure 
designed to maintain covariance. After a discussion of the general framework, we treat in detail 
the case of scalar fields, Yang-Mills gauge fields and the gravitational field. 



1. Introduction 

Despite the fact that the causal set program for quantum gravity pLj has been around for 
more than 20 years [21 El H]) identifying an appropriate formulation of the dynamics of the 
theory remains an open issue. Such a formulation would allow us to recover a smooth 
Lorentzian geometry with matter fields in a suitable continuum approximation, and would have 
a semiclassical limit in which general relativity coupled to those matter fields, including possible 
quantum corrections, is recovered. Some proposals have been made, most notably the classical 
sequential growth dynamics models of Rideout and Sorkin [5], a class of stochastic growth models 
defined by the fact that they satisfy causal-set versions of the general covariance and causality 
requirements. While a number of interesting results have been obtained on these models, no 
way to recover the continuum geometry and dynamics from them is known, and no proposal has 
been made for how to include matter fields in them. 

Here we will present a different approach to the dynamics of causal sets, one which is not 
based on general principles, but rather tries to mimick directly the continuum dynamics by 
using a discretized expression for the action, written in terms of variables that are meaningful 
for causal sets. Since a causal set is simply a locally finite partially ordered set, the only variables 
that characterize it are the number of elements in each region and the order in which they come. 
Thus, as far as gravity is concerned, our tasks are to express the Einstein-Hilbert action in terms 
of volumes of (finite) spacetime regions and causal relations, the continuum counterparts to the 
counting measure and the partial order, respectively, and then to discretize the result, making 
sure that we end up with a covariant expression. 



But one of our main motivations for proposing this approach is that the action for other fields 
can be treated in the same way, as we will show below for scalar and gauge fields. The reason 
for including matter is not simply completeness, or even a desire to develop the phenomenology 
of the theory, but the belief that ultimately the gravitational field will have to be treated on the 
same footing as the other forms of matter from the beginning. Since all matter fields will be 
coupled to the geometry at the causal-set level, and most causal sets are not well approximated 
by smooth manifolds, it may be that not only the metric structure of spacetime, but even the 
possibility of defining topological and differentiable structures depends on the matter fields, that 
their presence is necessary in order to understand how manifoldlike causal sets emerge. 

In section 2 we discuss the general framework and introduce the concept of Lagrangian 
generator, a building block from which suitable forms for the action of continuum fields can be 
written; section 3 treats the scalar field case and section 4 the case of a charged scalar field, i.e., 
one coupled to a gauge field; in section 5 we discuss the dynamics of the (Yang-Mills) gauge 
field itself, and in section 6 that of the gravitational field; the concluding section 7 contains 
some closing remarks. This article is based on the talk by LB at the DICE2008 conference, at 
which previous work along these lines |6] was presented. After the talk was given, however, RS 
developed an improved version of this approach, and we felt that it would be more useful for this 
writeup to reflect the recent developments. For this reason, the material in this paper reflects 
to a large extent the content of one chapter of the PhD dissertation by RS. 



2. Lagrangians, pre-Lagrangians and Lagrangian generators 

It turns out that the causal-set versions of the Lagrangians all look somewhat complicated, but 
they are all complicated in the same way. To emphasize the pattern, we introduce the concept 
of Lagrangian generator, together with a prescription for going from a Lagrangian generator to 
the corresponding Lagrangian. Our aim in doing that is to ensure that Lagrangian generators 
for all known flelds are simple-looking and the actual Lagrangians, as complicated as they might 
appear, can be "read off' from the expression for the Lagrangian generator. 

The key observation to guide us is the form in which one can write the scalars that appear 
in the usual continuum Lagrangian densities for the Klein-Gordon scalar, electromagnetic and 
gravitational fields. As will be shown in the next four sections, if a{p, q) := J~^{p) Ci J~{q) is the 
Alexandrov set defined by two causally-related points p ~< q m Minkowski space, then for any 
scalar fields and electromagnetic potentials varying linearly in inertial coordinates, respectively, 
the Lagrangian density at a point x £ a{p, q) can be written as 

d^cp d,^ = (^(g) - + 72^^;^ d'^- d-^^ - (2-1) 

/cem F'^'Ff,^ = / d'^r {a{p, r) + a(r, q) + a{q, p)) + 

vP; Q) J a{p,q) 

+ ^dw^ / + + ' (2-2) 

\Pi Q) J a{p,q) 

for appropriately adjusted pairs of coefficients {k,E), while for any metric whose expansion 
in Riemann normal coordinates only contains quadratic non-trivial terms (the lowest possible 
order, so we will still refer informally to this situation as the "linear case"), 

« = + ( „ lV(al^, r)) + V(a(r. ,))] ) . (2.3) 

for appropriately adjusted values for (/cgr, Egj.). Since a{p, q) is defined by causal relations alone 
without a direct reference to the continuum structure, and in particular observing that the right- 
hand sides of these integrals don't involve spacetime tensor indices, the replacement of integrals 
by sums makes it easy to generalize these expressions to a causal set. 



A point that needs to be emphasized is that the expressions ()2.ip - ()2.3p are covariantly defined 
scalars at a point, despite the fact that they involve the use of one Alexandrov set, as long as E 
is chosen to have the correct, dimension-dependent value. This, however, brings up a problem 
when we take these expressions over to the causal-set context: in light of the fact that the 
theory is formulated for arbitrary causal sets and not just manifoldlike ones, if E was viewed as 
a constant, it would be an extremely lucky coincidence that its value in the fundamental theory 
happened to be precisely the one needed for the exact cancellation of the non-covariant terms 
in the special case of a four-dimensional manifold. This, of course, is not satisfactory. 

For this reason, we will view E as a variable, subject to some physical laws that adjust 
its value to what it should be at any given point. The general idea is that the properties of 
a manifoldlike partial order determine the dimensionality of the manifold in which it can be 
faithfully embedded [Ij, which in turn determines the value of E. More specifically, while in 
the linear manifoldlike scenario the only field E depends on is gravity, as will be seen later in 
the non-linear case there is a slight dependence on the non-linear behavior of other fields, albeit 
very small. In the non-manifoldlike scenario, however, the dependence on non-gravitational fields 
might be arbitrarily large, just like the dependence on gravity would be. Consider a general 
theory for a set of fields of interest F ^ T, with a continuum Lagrangian density C{F, x, E) that 
can be obtained from expressions of the type (j2.ip ~ (j2.3p . Then 



is a function both of the point x and of the value of the fields F ^ T in its neighborhood. 
Notice that since E is not an independent field, a path integral for the field F will not include 
a separate integration over E. Rather, it will be of the form 



Let us now discuss how E{F, x) is determined. The definition of E for a general causal set 
requires a notion of degree of relativistic non-covariance for an expression such as (|2.ip - (|2.3p . 
which we will denote generically by C{F,E,p,q). If we identify a reference frame with the axis 
of an Alexandrov set (one defined by two points on the i-axis of the reference frame), then the 
degree of non-covariance represents the extent to which the value of C{F, E,p, q) depends on the 
choice of Alexandrov set a{p, q), and can be quantified as the range of variation of the values of 
C{F, E,p, q) as p and q are allowed to vary. 

In the case of linear fields, the non-covariance is zero as long as E is appropriately adjusted. 
In a curved Lorentzian geometry we can still have vanishing non-covariance, provided we go to 
the infinitesimal-a(p, q) limit. In the case of fields that can be seen as a discretized continuum 
it is not possible to select an infinitesimal region and the best one can do is constrain either the 
timelike length r(p, g) or the volume V{a{p,q)) between p and q to have a fixed (small) value, 
which means that some degree of non-covariance will appear, depending on the details of the 
discretization procedure; as long as the fields are linear or well-behaved functions that can be 
considered as slowly-varying, for an appropriately chosen E the non-covariance should be small. 

However, in the case of a general causal set, most of the assumptions that are made for 
the manifold can no longer be trusted, which means that it is possible for the degree of non- 
covariance to be large no matter what E is selected to be. However, it is still possible to formally 
select E in such a way that it minimizes the degree of non-covariance, even though the minimum 
might be very large. 

One feature of Lorentzian geometries that complicates this procedure is that if the manifold is 
not compact the set of pairs {p, q) may not be compact either, even when one imposes a condition 
on how far p and q are, so fields that are slowly varying in one Alexandrov set may appear to 







vary uncontrollably in another, highly boosted one that occupies a stretched region "near the 
light cone of point p". Therefore, if all possible Alexandrov sets a(j),q) with T{p,q) = tq were 
considered, the value of E might have nothing to do with the one it would have had in the linear 
case. For that reason, we constrain the varying Alexandrov set to vary within the boundaries of 
some other, larger "fiducial" Alexandrov set a{P,Q), with 

n=Tip,q)<TiP,Q)=T2. (2.6) 

The variation, then, is defined by 

YSiir,{^,F,E,P,Q) ■.= ma^{C{F,E,p,q)\ P ^p^q^Q A T{p,q)=n}, (2.7) 

and for any point x in a causal set, ari,T2 (-^; is defined as the set of triples {E, P, Q) for which 
var,-^ (£, F, E, P, Q) is minimized with the constraint that t{P, Q) = T2: 

ar,,r2iF, x) = {{E, P, Q) I r(P, Q) = T2 A V E' , P' , Q' with t{P', Q') = T2, 

Ya.Tr, {C, F, E', P', Q') > vavr, (£, F, E, P, Q)} . (2.8) 

Typically, aT-^^T-2(-^) ^) is a one-element set, and the Lagrangian at x can simply be taken to 
be C{F, E, P,Q), where {E,P,Q) is the unique element of that set. But in order to formally 
accommodate the cases where ari,r2(-^) ^) has more than one element, the Lagrangian at x is 
formally defined as the average of the above over all elements of aT-i,T2(-^) a;): 

Cr,,rAF,x)= \ Yl mE,P,Q), (2.9) 

where jj{ } denotes the cardinality of a set. 

Now, dropping the issue of the value of E and the selection of p and q for a second, let's switch 
gears and go back to the expression for C{F, E,p, q), which, from now on, will be referred to as 
a "pre-Lagrangian" . As the examples in the beginning of this section illustrate, the expressions 
for C{F,E,p,q) for scalar, gauge, and gravitational fields look similar, in the sense that they 
are all linear combinations of a function of p and q and some form of integral of the same 
function over the interior of a{p,q)- So the natural question arises: why do Lagrangians take 
this particular form and not some other one? That question is answered by introducing the 
concept of Lagrangian generator and formally defining a procedure of going from a Lagrangian 
generator to a pre-Lagrangian in such a way that the above mentioned linear combination arises 
in a natural way if one formally follows the steps of the procedure. 

For simplicity, let's consider first an expression for a pre-Lagrangian of the general type ()2.3p . 



Cj{F,E,p,q)=J{p,q) + -j^ I d^rj(p,r), (2.10) 

where ^7 is a real- valued function depending on the set of fields F, which maps {p, q) 
JiyF^p, q) G M. After the selection of E and q will be made, the value of the pre-Lagrangian for 
that specific selection will be identified with the actual Lagrangian density at p. To put (j2.10p 
in a useful form, define two maps f,g:S^^S'^, where here S is any set, as follows: 

f{a,b,c) = {a,b), g{a,b,c) = {a,c) . (2.11) 

With these definitions, 

[ d''rJ{FJ{p,q,r))= [ d'^r J{F,p,q) = V{a{p,q)) J{F,p,q) (2.12) 

Ja(j),q) Ja{p,q) 

and 



[ d''rJiF,gip,q,r))= [ d^rJ{F,p,r), (2.13) 

Ja(p,n) -JaipM) 



i(p,ij) Ja{p,g) 
and a pre-Lagrangian equivalent to (j2.10p in the continuum is given by 

£j{F,E,p,q)= ] [ d''r[j{F,f{p,q,r))+Ej{F,g{p,q,r))], (2.14) 

V [Ol\p, q) ) J a{p,q) 

because V{a{p,q)) and r*^ are proportional, through a constant kd (see (j3.6|) below) that can be 
absorbed into E. In fact, since T{p,q) = ti will also be a constant, V{a{p,q)) may be omitted, 
and in order to preserve time-reversal symmetry, we would like to replace J^iF, f(p,q,r)) 
and J{F,g{p,q,r)) by J{F, /{p, q,r)) + J^F, f{q,p,r)) and J{F, g{p, q,r)) + J{F, giq,p,r)), 
respectively, up to an unimportant factor of 2. Thus, the pre-Lagrangian can be written as 



£j{F,E,p,q)= [ d'r (j{F,fip,q,r)) + j{F,f{q,p,r)) + 

+ E[j{F,g{p,q,r)) +j{F,g{q,p,r))]) . (2.15) 



The Alexandrov-set function J'(F,p,q) will be called a "Lagrangian generator" for the pre- 
Lagrangian p.lSp . It is the basic object in the dynamics of the field F, and needs to be defined 
without reference to spacetime tensors or a differentiable structure on the underlying set. 

More generally, in order to allow Cj to include multiple integrals such as in ()2.ip and ()2.2p , 
f{p, q, r) and g{p, q, r) should be generalized to f{p, q, ri, r„) and g{p, q, ri, ...,rn)- For the sake 
of completeness and convenience to the reader, our final definition of Lagrangian generator J' 
will include both the prescription for going from the Lagrangian generator to the pre-Lagrangian 
JC-j{F, E,p, q) just discussed, and the transition from the pre-Lagrangian to the actual pointwise 
Lagrangian Cj{F,x), discussed earlier in this section: 

Definition: Let J- be the set of possible distributions of a set of field of interest. A Lagrangian 
generator is a triple {J, /, g) with J : J" x 5"^ ^ M, and /, g : S^+r" ^ s'^, with n<2 + m. The 
pre-Lagrangian corresponding to this Lagrangian generator is Cj : J- x x S'^ ^ given by 



p^Ti^q 



in the causal set case, with the sums replaced by integrals in the continuum case. The "variation" 
of this pre-Lagrangian in a(P, Q) is the real- valued function defined by 

var,, {J, F, E, P, Q) = max{£^(F, E,p,q) \ P < p < q < Q A t{p, q) = n} , (2.17) 

and for any x G S, the (F, ri, r2)-based fiducial neighborhood of x is given by 

ar,,r^{F,x) = {{E,P,Q) I t(P,Q) = ra A yE',P',Q' with t{P',Q') = ra, 

vavr, (£, F, E', P', Q') > vairi {C, F, E, P, Q)} . (2.18) 

Finally, the pointwise Lagrangian density corresponding to is given by 

' E C(F,E,P,Q). (2.19) 

HI ri,T2V ; E,P,Qea{x) 

Let us now see what this expression becomes for some important types of fields. 



Cj{F,E,p,q)= ^ i^j{FJ{p,q,n,...,r„^))+j{FJ{q,p,n,...,rm)) + 

+ E[j{F,g{p,q,n,...,rm)) +j{F,g{q,p,n,...,rm))]) , (2.16) 



3. Scalar Fields 

The Lagrangian generator for a scalar field (f) is given by {J',f,g), where 

Ji^, r, s) = {<P{r) - - i mV'(r) (3.1) 

/(n,r'2,r3,r4) = (ri,r2) , 9{n,r2,r3,u) = {r3,r4) ■ (3.2) 

For reasons that will soon become apparent, m is not the actual mass, although it is related to 
it. In fact, m is assumed to be very small, of order r. The above expression implies that the 
pre-Lagrangian for a scalar field is given by 

jCict),E,p,q)= J d''rd''s{j{f{p,q,r,s))+EJig{p,q,r,s))) (3.3) 

= J d\d'^s{j{cP,p,q) + EJ{cP,r,s)) = J{cP,p,q)V\a{p,q)) + E J d\ d'^s Ji<P,r, s) . 

By remembering that the volume of an n-dimensional ball is 



2„(d-l)/2 i't/2 d-l 

V{a{p,q)) = ,,,,, / ( - - |t| ) dt = kaT\ (3.5) 



we obtain, for the volume of an Alexandrov set in Minkowski space, 

(d-i)r((d-i)/2)y_,/2V2 

where 

27r('^-^)/2 

^'^ " (d-i)r(2'^-2d(d-i)) ■ ^^-^^ 

Substituting the above expression for volume, along with the Lagrangian generator J', in the 
expression for C, we obtain 

C{ci>,E,p,q) = kjT^''ip,q)[{^{q) - ^{p))^ -\m^ {cl>\p) + 4>\q))] 

+ E J d'^rd'^s [((^(r) - (l){s)f - m'^^\r)] . (3.7) 

We would now like to compute the Lagrangian density, if the spacetime is assumed to be flat 

Minkowskian and (f) is assumed to be linear. 
Let's start from the mass term. It is given by 

£m{<P,E,p,q) = !I^r'"'ip,q) [cf" {p) + <t>\q)) + E j d'^r d'^s ^\t) . (3.8) 

The above expression tells us that the leading order of the mass term is m? r^'^. When we will 
get to the kinetic term, it will be shown that the order of magnitude that we are interested in 
is r^'^^^. Thus, if m is assumed to be of the order of r, then the leading order is r^'^^^ which 
coincides with the leading order for the kinetic term. 

This means that, as far as the mass term is concerned, we can throw away all the higher-order 
terms. This can be done by using the approximation ^ k. ^q, which tells us that 

Cm{^, E,p) ^ rr? kj t''%, q) 4'1 + Em'' cf>l V\a{p, q)) 

= {l + E)mHj<f>lT^''{p,q). (3.9) 



Now let's look at the kinetic term. The linearity assumptions imply that in the integrand 

(0(r) - cj){s)Y = K - sn {r" - s'')d^(j)d,4>. (3.10) 

Consider a coordinate system in which the t axis passes through p and g, while the origin lies 
midway between these points. Denoting t(p, q) by r, 

p= (-r/2, 0,0,0), g= (r/2, 0,0,0). (3.11) 

In this coordinate system, 

[ d'^rr'^ = 0, (3.12) 

Ja{p,q) 

since the above integrand is antisymmetric with respect to the center of the Alexandrov set. By 
slicing the Alexandrov set into t = constant balls, we get 

where 

27r('^-i)/2 

" {d - 1) r(2'^-i {d-l)d{d + 1) {d + 2)) ■ ^^'"^^^ 
Furthermore, it can be shown that 

(x'=)2d'^^ = /rflr'^+^ (3.15) 



/ 

Ja( 



la{p,q) 

where, by cylindrical symmetry, the coefficient is the same for each k and is given as 

27r<i/2-i 

= (d-2)r(d/2-i) • ^^-^^^ 

Substituting these expressions into the integral in i2kin we obtain 

/„ d—l „ d—1 

d'^r d'^s {<P{r) - <P{s)Y = n d'^s^ d^r {r'^f +U d\^ ^ d'^s {sf'f 

= 2 kd t'"'+'' [ [ho + hi id - 1)) {doct>f - hi {d - 1) d^'<j) d^,<j)] . (3.17) 

Thus, the kinetic term of the pre-Lagrangian is 

jCun{^,E,p,q) = kdT^^+^ {{do^f [kd+2Ed {Ido+hi {d-l))]-2Edhi {d-1) d^cl>d^<p} . (3.18) 

Switching from the coordinate system in which the t axis passes through p and q to an arbitrary 
one, the result becomes 

£kin(<^, E, p, q) = kd T^'^ip, q) {q^ - p^) {q" - p"") d^cj) 8,^ [ka + 2 Ea {ho + hi {d - 1)) ] + 

- 2kdEdhiT^''^^{d-l)d^'(t>d^(t>. (3.19) 

If the choice of points p and q varies with the constraints that both the midpoint between 
p and q as well as the Lorentzian distance between the two points are fixed, then £{(/), E,p,q) 
undergoes a variation of order r^*^"*"^ due to the (g'^ — p^^) {q" — p") 8^(1)8,^(1) term. The mass 



term, on the other hand, only gives variations to higher orders. Thus, if variations of orders 
higher than r^'^^^ are neglected, then the variation can be "minimized", or in this case set to 
0, if the first term in (13.190 (depending on the "tilting" of the axis of the Alexandrov set, while 
the second one is constant) vanishes, or 

kd + 2Ed{ldo + Idi{d-l)) =0, (3.20) 

which determines the value of Ed'. 

^'^ " ~ 2 {Ido + Idiid - 1)) ■ ^^"^^^ 
Substituting this into the expression for the Lagrangian gives 

IdQ + Idi{d-l) ^ d y 2{ldo + Idi{d-l))J ^ ' 

This Lagrangian can be rewritten as 

C = h'/-^d'^^dd,^d-'^cPl (3.23) 
where vo is the volume taken up by one point and 



hd^d = kr'^+^ J^I^li^^^ ^ (3.24) 
^^0 V Ido + /di (d - 1) ^ ^ ^ 



md = ^T\ 1-— -m. (3.25) 

vo y 2{ldo + Idi{d-l)) 

In future sections other fields will be similarly scaled, but the coefficients will be different from 
field to field. At first this might seem wrong since the kinetic terms of all Lagrangians have 
the same coefficient 1 in standard quantum field theory. But it can be easily shown that this 
difference does not amount to anything but a change of overall factor: 

J [P(/>i...P</>n] exp{iS{<Pi) + ... + iS{<Pn)} (3.26) 

= Pld ■■■Pnd / [V(l)id...'D(l)nd] exp {\S{(f)id/p) + ... + lS{<i)nd/ P)] , 



where 

(pkd = Pkd(pk- (3.27) 



4. Charged Scalar Field 

Consider a charged spin-0 particle, described by a set of complex scalar fields (j) = {'Pit--! 4'n) 
coupled to a SU(n) gauge field. (Notice that in this approach to the dynamics of matter fields in 
causal set theory, although it will be assumed that spacetime is discretized, the internal degrees 
of freedom will still have a continuous invariance group.) In the continuum, the dynamics of 
such a field can be described starting with the matter Lagrangian density 

Cr^ig^,, <P, A^; x) = \ \g\^'^ [g^'' {D^4>)^ {D,<^) - , (4.1) 



where the gauge-covariant derivative is defined as usual by D^cf)"" := d^(f>°' + ie Afj,°'i, cj)^ , and 
= A^j^ is the Lie-algebra- valued connection form representing the gauge field on a 
differentiable manifold. (Here, Latin indices a, h, are Lie-algebra tensor indices, while /c, 
I, label elements of the basis of the Lie algebra.) In the causal set context, the scalar field 
will be simply replaced by a corresponding field defined at each causal set element, but to write 
down the action it is important to specify what variables will replace A^. 

A gauge field is defined in terms of holonomies, where the word "holonomy" refers to the 
group transformation corresponding to the parallel transport of a Lie-algebra-valued field such 
as (j)"" between two points p and In a differentiable Lorentzian manifold M (of dimension d), 
holonomy is defined as a function a : M x M — > TSU(n), where TSU(n) is the tangent bundle 
to SU(n), and consists therefore of all traceless n x n tensors. This map assigns to any two 
elements p, q € M the holonomy of along the geodesic segment 7(p, q) connecting p and q 
in M, given by the path-ordered exponential 

a{p,q)= [ V^'dx^ (4.2) 

in terms of which the expression Dfj_(j){x) appearing in the scalar field Lagrangian arises from 
the leading-order term in the expansion of the expression (1 -|- a{x,y)) {4>{y) — <P{x)) ■ 

This means that the causal set version of the charged scalar field Lagrangian can be obtained 
by making some simple substitutions in the one obtained in the previous section for the Klein- 
Gordon field. Thus, if the gauge field a is assumed to be fixed and not subject to any Lagrangians, 
then the Lagrangian generator for a matter field G interacting with a is given by (J7, /, 5), 
where 

Jsc(0, a, r, s) = W) - a(r, s) (/.(s)p + \ (</.* (r) + 0*(s)) (0(r) + (/>(s)) (4.3) 

and 

/sc(n,'^2,?'3) = (n,'^3) , 9sc{ri,r2,n) = {ri,r2) . (4.4) 
However, as discussed in the next section, a itself is subject to a Lagrangian generator given by 

JYM{a,ri,r2,r^) = tr[(a(ri,r2) + a(r2,r3) + a(r3,ri))^] , 

9YM{ri,r2,r^,ri,r^) = {ri,r2,r^) . (4.5) 

In order for the theory to possess an SU(n) symmetry in the setup, these two Lagrangian 
generators are combined into one as (i7tot , f,g)i 

Jtot(a,ri,r2,r3) = Jsc{4>,a,ri,r2) + JYM(a, n, '^2, "rs) , 
/tot(?'l,?^2,?'3,'^4,?^5) = /YM(n,r2,r3,r4,r5) , 

gtot{ri,r2,r3,r4,r5) = 5YM(n, ^2, r3, r4, rs) . (4.6) 

5. Yang-Mills field 

In this section, the main goal is to express the Yang-Mills Lagrangian density, 

CYMig^..,A^■,x) = l\g\^/Hr{F^,F'"'), (5.1) 

in terms of the holonomies for the gauge field introduced in the previous section, as well as 
variables describing the geometry that are meaningful for causal set [7j. Once this is done, the 
Lagrangian density can be easily rewritten in the discrete setting. 



In the causal set context, a gauge field is defined as a map a: ^ TSU(n). The Lagrangian 
generator for the gauge field is given by {J,f,g), where 

J{a, n, r2, rs) = tr[(a(ri, + a(r2, rs) + a(r3, n))^] (5.2) 

and 

f{ri,r2,r3,r4,rr,) = (ri,r2,r3) , c/(ri, r2, rs, r4, rs) = (r3,r4,r5) . (5.3) 
This means that in the case of Minkowski space the pre-Lagrangian is given by 

£(o, E, p,q)= I d'^r d'^s [j{a, f{p, q, r, s,t)) + J {a, /{q, p, r, s,t)) + 

+ E {J{a,g{p,q,r,s,t)) + J{a, g{q,p,r, s,t))] 
= I d'^rd'^sdh[j[a,p,q,r) + J{a,q,p,r)+E{J(a,r,s,t)+J[a,r,s,t))\. (5.4) 

Jr(p,(j) 

Integrating (trivially) over s and t in the first two terms, and using the fact that J{a,p,q,r) = 
^{a, q,p,r), we get that 

C{a,E,p,q) = 2(v^{a{p,q)) f d'^r J{a,p,q,r) + E f d'^r d'^s dH J{a,r, s,t)y (5.5) 

JT{p,q) ■Jr{p,q) 

Assume that the gauge field is differentiable and reasonably well behaved. In particular, it is 
well-behaved-enough for -F^i/^ to be approximately constant in a{P, Q) whenever {E, P, Q) G 
o^ri,T2('^) ^^) for some u E 5. Assume for definiteness that the three points are spacelike related. 
Choose a coordinate system so that r coincides with the origin, the x axis points from r to s, 
and the y axis is perpendicular to the x axis in the r-s-t plane. Then in this coordinate system 
a = (0,0,0,...), 6= (0,61,0,...), c= (0,ci,c2,...). 



The flux of Fnu'' through the interior of that triangle is expressed by the relationship 



a(r, s) + a(s, t) + a{t, r) = ^ Fi2^ + ... (5.6) 
This result generalizes to points at arbitrary locations, and can be written covariantly as 

a(r, s) + a{s, t) + a{t, r) = \ F^,'' T>' (s^ - r^^) - r'') + ... (5.7) 
Recalling that, for SU(n), tr(r'^T') = C2 5fc^ to leading order in the separation between points, 

tr[(a(r, s) + a{s, t) + a{t, T)f] = ^ F^,'' (s'^ - (f - r^) F^/ - a") {c" - a'') . (5.8) 
Let's start from 

/ d'^rd'^sd'^tj{a,r,s,t). (5.9) 

Expand the right-hand side of Eq (j5.8|) . and integrate term by term. Clearly any term with 
an odd number of powers of any variable will integrate to 0. Thus, the only terms that may 
potentially survive the integration are those of the form r^ r'^ r" or quadratic terms in two of 
the three points. Simple counting of terms gives 

I d'^r d'^s d'^t tr [(a(r, s) + a{s, t) + a{t, r)f] 

Jp~ir,s,t^q 

= ^F^^^Fp^^ j d'^ad%d'^c{s>' -r>'){t^ -r''){sP -rP){t'' -r'') 
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a(r^)2)( / d%is''f] + 

a{p,q) 



V^F^^^Fp/ / d'^ar'^r'^r^'r'^ 

Ja(p,q) 



(5.10) 



where V is the volume of the Alexandrov set a{p, q). 

The only terms of F^J^ Fp^^ r'^ that survive integration are the ones whose indices 

are pairwise equal. But if either fi = u or p = a, then F^i/'^ = or Fpcr^ = 0, respectively, which 
would set the whole thing to 0. Thus, the only options are = p, v = a and u = p, p = a. The 
antisymmetry of Fp,y^ then implies that these two cases are opposites of each other, which in 
turn imphes that F^^'' Fp„^ r^" r" rP = 0. Thus, Eq ([KTO]) b ecomes 

f d'^r d'^s dH tr [(a(r, s) + a(s, t) + a{t, r)f] 



Q 7 3d+4 d—1 

^ '^l" E {J'J' Ei^^o'f + {j'r , (5.11) 



k i=l i<j 

where JP = r^'^^^ f^^^ d'^x (x'^)^, or in other words 

o_(d-l)/2 

Now let's move to the second integral, ja{pq) d'^a^tr[(a(p, x) + a{x,q) + a{q,p))'^], where p < q 
are the endpoints of the Alexandrov set. Rewriting Eq (j5.8p in terms of the points p, x, and g 
gives 

tr [(a(p, x) + o(x, q) + /(g, p))^] (5.14) 

= \ C2 Fp/ {p^ - x^) (g- - x'^) Fp^'^ {pP - xP) {q- - x^ . 

Again this can be expanded and integrated term by term. There are several conditions each 
term has to meet, in order for its integral not to vanish. First of all, it needs to contain an even 
number of factors of x. Secondly, as was shown before, for symmetry reasons 

Fpy^Fp^^ I d'^x xP x" xP x" = . (5.15) 

Finally, Fp^pPp" = Fp^qPq" = 0, and the identities p = (-§,0,0,0) and q = (§, 0,0,0) imply 
FpypPq'^ = —FpypPp^ = 0. The only terms in Eq ()5.15p that do not vanish for any of the above 
reasons are 

Fpu"" Fp/ pP x^ pP x-^ , F^/ Fp/ pP x^ xP q'' , 
Fpu^ Fp„^ xP pP X- , FpJ' FpJ' xP xP q'' . 

Plugging in the coordinate values of p and each of the above four expressions evaluates to 
\ t"^ Fpo' Fpo^ xP xP. In order for this not to be an odd function, p = p has to hold, and in order 
for FpQ to be non-zero /U / has to hold. Thus, this becomes | {Fio^y (x*)^ and, since there 
are four such terms, the integral becomes 

d'^x tr[(a(p, x) + a{x,q) + a(g,p))^] 

a{p,q) 

i C2 E / d'^x {F.o'f {xr = I C2 r^+'j' Y^i^io"? , (5.16) 



where, based on rotational symmetry, = ... = J'^~^ was used. 

Substituting this into the original expression for the pre-Lagrangian, and doing some basic 
algebra, the latter becomes 

C{a,E,p,q)=T^''+^C2kdJ^ [iSEj' + lkd)^iFtof + 3Ej'Yli^tjf) ■ (5-17) 

i<j 

In order to get rid of the variation that results from different choices of Alexandrov sets, we 
would like the above expression to be relativistically invariant; in other words, we would like it 
to be proportional to F^'^F^i,. This means 

3EaJ'' + ^kd = -3EaJ\ (5.18) 

which implies 

Substituting this expression for E^ gives the Lagrangian 
which can be expressed as 

£(a,p,g) = -^F/-Fd^,, (5.21) 
where vq is the volume taken up by a single point if F^'^'^ is defined as 



F.- = ».7'F-J-^1^. (5.22) 



This means that 



From this we can determine how the charge changes from dimension to dimension. At a first 
glance, since we haven't yet defined a Lagrangian for fcrmions, we are only ready to talk about 
the charge of bosonic fields. However, a simple symmetry consideration allows us to overcome 
this barrier and include fermionic charges in the discussion. Whether a field is bosonic or 
fermionic, we would like to be able to say 

d^^d^, + eA^. (5.24) 

We would also like to be able to say 

^^,^^^ + edAd^. (5.25) 

This means that, regardless whether the field in question is bosonic or fermionic, and regardless 
of any other properties of the field (such as mass) it should satisfy 

eA^ = edA^. (5.26) 

This immediately implies that, both for bosons and fermions, 



6. Gravitational field 

The Lagrangian generator for gravity is the triple ( JT", /, g) , defined simply by 

^, , [(SttG)"^ iir^t^s _ , , , , , , , 

J{r,s,t) = < . , f{r,s,t,u) = {r,s,t) , g{r, s,t,u) = {r,t,u) . 

I otherwise 

(6.1) 

This means that the pre-Lagrangian is given by 

£ = J d^r d'^s VI det g{r) | VI det g{s) \ [j {f{p, q, r, s)) + E j{g{p,q,r,s))] . (6.2) 

Since J" is a constant inside a certain domain and is outside, the above integrals amount to 
restricting r and s to certain domains. The first term is J'{p, q, r), thus its restriction is p ^ r ^ q. 
That statement, of course, is trivial, which means that the first term has no restriction at all. 
On the other hand, the second term is J^{p, r, s), thus it has a restriction p -< r -< s -< q. Taking 
this into account, the pre-Lagrangian becomes 

^ = ^( [ d''rd''s^/\detg{r)\^/\detg{s)\+E [ d'' r d'' s ^/\ det g{r)\ ^/\ det g{s)\) 

= {V{a{p,q))f + E f d''s^\detg{s)\V{a{p,s)). (6.3) 

Ja{p,q) 

One can parametrize the interior of the Alexandrov set a{p, q) with normal geodesic coordinates 
around p. Suppose that Greek indices /x, u are used for coordinates satisfying 

9i^i^{p) = ViJ-u , (6-4) 

but which are otherwise arbitrary, not necessarily geodesic. Then one can define the normal 
geodesic coordinates derived from the above in the following way: 

1/2 



r'^ = r,/^-5,7prij,(/ SM.dx'^dx'^) , (6.5) 



where 'jab denotes the geodesic segment connecting a and b. Prom this, it is straightforward to 
see that the following equation is satisfied exactly: 

r(p, r) = ijfipr'^r" . (6.6) 

It should be noticed, however, that the above is true only if one of the two points is p, while 

r{r,s)y^ViiAs^-rn{s'-rn. (6.7) 

However, a different causal relation, -<p, will be introduced in addition to the already existing 
one, -<. While -< matches flat spacetime expectations only if the pair of points of interest includes 
p, -<p does so for arbitrary pairs of points: 

a^ph ^ riup {If - a^) {}f - a^) > . (6.8) 

At the same time, we will retain the original causal relation, -<, for which the above is not true: 

^ 7^ . (6.9) 



a{p,q) 



Let's introduce the following notation: 

J+(a) = {r y a} , Jp{a) = {r >-p a} , 

J- (a) = {r -< a} , Jp{a) = {r -<p a} (6.10) 

a{a,b) = J+{a)nJ~{b), ap{a,b) = J+{a) n J~{b) , (6.11) 

and define Vp{a,b), not to be confused with V{ap{a,b)), as follows 

Vp{a, b) = ka (?7m. {h^ - an (b" - a^f' . (6.12) 
It is important to notice that 

voliapia,b))^Vpia,b), (6.13) 

because the definition of Vp neglects the l^l^^^ factor in the volume element. But, for our 
purposes, Vp as defined above is the simplest one to use. On the other hand, V{a, b) is defined 
as the actual volume of the Alexandrov set, 

V{a,b) =vol{a{a,b)). (6.14) 

Then, the integral of the volume can be expanded as follows 

d'^r\/|det5(r)| V{a{p, r)) 
= [ d''rVp{p,q) + Auip,q)+A2d{p,q) + A3a{p,q) + 0{T^''+n, (6-15) 

Jap(p,q) 

where Aifi{p,q) is the correction coming from the mismatch between a{p,q) and ap{p,q), 
A2dip,q) is the correction due to the error in V{a{p,r)), and A3 is the correction due to the 
use of d^r instead of d'^r^^ (— 1)"^"^ det g for the volume element. These corrections are formally 
defined as follows: 

Aid{p,q)= [ d<'rVp{p,r)- [ d''rVp{p,r) (6.16) 

A2d(p, q)= [ d^r {V{p, r) - Vp{p, r)) (6.17) 
^3dip,q)= [ d^r{y/i-iy-^deig-l). (6.18) 

Whenever any of these three correction terms is computed the other two are neglected, since 
the "correction of the correction"' is of order r^^^^^, while the calculation is performed to order 
^2d+2_ 'j']2us^ th e shape of the Alexandrov set is assumed to be unchanged in the calculation of 
A2d and A^^, ^/ (— 1)<^~^ det 3 is dropped in the calculation of Ai^ and A2d, and the correction 
to V{a{p, x)) is neglected in the calculation of Ai,;; and A^d- 

Let's start by computing Ai^. In normal coordinates, the light cone of p is not deformed by 
curvature, while the light cone of q still is, 

J^ip) = 4{P) , J'iq) + Jp{q) ■ (6.19) 
Substitution of above into simple set theory algebra gives 

ap(jp, q) \ aip, q) = J+ (jp) n (J- {q) \ J' (q)) C J' (q) \ J" (q) (6.20) 

aip, q) \ apip, q) = J+{p) D (J-(g) \ J~iq)) C J'iq) \ J;{q) . (6.21) 



Assuming that Jp{q) and J~iq) are very close to each other, the above imphes that most of the 
contribution to Ai,^ comes from the vicinity of Jp{q)- So, let's evaluate Vp{p,r) for r G Jp{q)- 
In other words, assume that 

This implies that 

r^^^ = {r^ - (r - r^f = r (2r° - r) , 

which, using the fact that 



(6.22) 
(6.23) 



2r" - r = r" - (r - r°) = r° - JEi^'^V , 



becomes 



Now define two functions xi and X2 as follows: 



Xiir) 
X2ir) = / 

J a 



a{r,q) 



To 



d'^s kd iVfip s 

d/2 



H „u\d/2 



(6.24) 
(6.25) 

(6.26) 



d'^skd 



To 



- . (6.27) 



As a consequence of the fact that J'^{p) = Jp{p), we get 

Xi{p) = X2{p), 

and it is easy to see that 



XiiQ) =X2iq) = 0. 



Therefore 



^^''\^dx2{riT)) 



nd = Xiip) = -iX2iq) - X2ip)) = - ~. 
where r(r) is defined as a point on jpg whose distance to r is r: 

7pr(T) C 7pg , l{lpr{T)) = T ■ 

This immediately implies that 

St > 0^ r{T) ~< r{T + 5t) . 
Furthermore, in normal coordinates 



(6.28) 
(6.29) 
(6.30) 

(6.31) 

(6.32) 
(6.33) 



where V^^ is a tangent vector to ^pq at p. The fact that the latter is timelike implies that 

5t>Q ^ r(r) ^p r{T + 5t) . (6.34) 

In order to compute 6.X2{t{t))/(\t notice that, schematically, 

SX2{r)= [ f- [ f- [ f+ [ f= [ f-[ f-j f+[ f 



{S{\S2)\{S^\Sa) 



{S-i\Si)\{S^\S2) 



(S3\S4)\(S4\S3) 



{Si\S3)\{S'j,\Si) 



/, (6.35) 



where 

Si = a{r{T + 5T),q) , S2 = ap{r{T + 5T),q) , S3 = a{r{T), q) , S^ = ap{r{T),q) . (6.36) 

Substituting both r(r) -< r(r + 6t) and r(r) r{T + 6t) and using simple set-theory algebra, 
one obtains 

a{r{T), q) \ ap{r{T),q) \ (a(r(r + (5r), g) \ ap{r{T + (5r), g)) 

= (J+(r(r)) \ J+(r(r + <5r))) n {J-{q) \ ap(r(r, <?))) C J+(r(r)) \ J+(r(r + <5r)) (6.37) 
(^pirir), q) \ a{r{T),q) \ {ap{r{T + 6t), q) \ a{r{T + (5r), g)) 

= (J+(r(r)) \ Jp+(r(r + 5r))) n (J-(g) \ a(r(r, g))) C j;(r(r)) \ J+(r(r + 6t)) (6.38) 
(a(r(r + (5t), q) \ ap{r{T + (5t), q)) \ (a(r(r), g) \ ap(r(r), g)) 

= a(r(r + 6T),q) n (Jp+(r(r)) \ J+(r(r + 6t))) C J+(r(r)) \ j;(r(r + dr)) (6.39) 
(ap(r(r + 6t), q) \ a{r{T + 6t), q)) \ {ap{r{T), q) \ a(r(r), g)) 

= ap(r(r + 6t), q) n (J+(r(r)) \ J+(r(r + 5t))) C J+(r(r)) \ J+(r(r + dr)) . (6.40) 

Thus, all four integrals are performed either over a subset of J"^(r(r)) \ J+(r(r + 5t)) or over a 
subset of Jp{r{T)) \ Jp{r{T + St)). In either case, the range of integration is in the vicinity of 
the lightcone of r(r). 

Now let s be an arbitrary point in that region. Remembering that r denotes the distance 
from p to the point of intersection of that region with the geodesic 7pg, and is not to be confused 
with the distance from p to a floating point, we have 



s 



e [J+(r(r)) \ J+(r(r + 6t))] U [j;(r(r)) \ j;(r(r + 6t))] 



d/2 



/J2{i^)2 = - r + 0(r2) ^ k, [to [s^ - = {tt^)^'^ + OiT^^') . (6.41) 

Thus, all integrals reduce to integrations over a constant: 
/ fd'^x- I fd'^x- I fd'^x+ 1 fd'^x 



-/ -/ 



+ / fd^'x 

{Si\S2)\{S3\S4) J {S3\S4)\{Si\S2) J iS3\S4)\iS4\S3) J {Si\S3)\{S3\Si) 

irro)'/'([ -f -[ +[ )d'x 

^ J {Sl\S2)\iS3\Si) J {S3\S4)\iSi\S2) J {S3\S4,)\{S4\S3) J {S4\S3\{S3\S4) ^ 



(rro)'^^^ if d'^x- f d'^x- f d'^x + f d'^ 
^Js, JS9 J S3 J Si 



'Si J S2 J S3 J S4 

= (tto)"/^ {V{Si) - V{S2) - V{Ss) + ViS^)) . (6.42) 

Now, V{Si), V{S2), ^('S'3) and Vi^S^) can be read off from reference [8j as follows: 

V{Si) = [kd + {Ad Rg^, + Bd R^,) {q^ - (r + 5t)) {q"" - r"" {t + 6t))] T'^(r(r + 6t), q) 
V{S2) = kdT''{r{T + 6T),q) 

ViSs) = {kd + {A Rg,, + B R,,) {q^ - t^{t)) {q'^ - r'^{T))) T''{r{T),q) 

V{S4) = kdT''{r{T),q). (6.43) 



Since the coordinate system is defined in terms of geodesies coming out oip, all of these geodesies, 
including jpq are, by definition, straight lines in the chosen coordinate system. Therefore, it can 
be assumed that coincides with the t axis, which simplifies the above equations: 



V{Si) =kd{Ti-T- drf + {AdR + Bd Roo) {n - t - Srf^^ (6.44) 

^(^2) =kd{n-T- 6Tf (6.45) 

V{Ss) = kd (ri - t)^ + iAdR + Bd Roo) in - r)'^+2 (6.46) 

V{S4) = kd (ri - t)^ . (6.47) 

This implies that 

X2(r(r + St)) - X2(r(r)) = {nT)''/^ViSi) - ¥(82) - ViS^) + ViS^)) (6.48) 

= {AdR + Bd i?oo) (riT)'^/2 d_ _ ^^d+2 ^ Q(^d+2 ^^^)2) ^ 

which in turns implies that 

^^^^ = {AdR + Bd Roo) {nr)'/' ^ (n - r)'^^^ . (6.49) 

Thus, 

Aid = Xi{p)-Xi{q) = X2{p)-X2{q) = -{AdR + BdR^) £ dr (rir)''/^ i.(^^_^)<i+2 . (6.50) 
The binomial expansion of (ti — r)'^+^ gives 

Ai. = -{Ad R + Bd Roo) £ dr {r,r)'/' A |^ + 2^J _ ^^^^^ 



Evaluating the derivative and combining it with the t'^/'^ factor we find that 

Ai. = -{AdR + Bd Roo) £ Arrt" ( Y.^-!)' X ^) ^1'^'"' ^ r^"^'-') ■ (6.52) 



■A;=0 

Finally, integration of this last expression gives 

d+2 



d + 2\ k 



A,d ={AdR + Bd Roo) r^' ("^ 1 ^) 



d/2 + k 



(6.53) 



Now let's compute ^24- We will go back to our usual notation, where r (which will be denoted 
by X in order to avoid conflict with r being a radius of a ball) is an arbitrary element of a{p, q), 
and it is no longer assumed to lie on 7^^. As was stated earlier, the A2d correction is due to 
the error of computing the volume of a{p, x), where x is far away from the boundary of a{p, q). 
Since that correction is already of order r'^"'"^, its integral over a{p, q) is of order r^'^"'"^, while its 
integral over the "corrections" to the shape of a{p, q) is of order r^'^"'"^. For this reason, the latter 
term will be neglected, and it will be assumed that the shape of a{p, q) has not been affected 
by curvature, while the shapes of a{p,x) has. Furthermore, in light of the geodesic coordinates. 



while the corrections to the volume of a{p, x) are not neglected, the flat space equation for the 
distance will be used in computing them. Thus, the equation for becomes 

By using the cylindrical symmetry, this becomes 

^2d = Ad 

f d'^xix^fix^Xpf/'^, 

where, due to cylindrical symmetry, (x'^)^ was replaced with (x^)'^. 
Prom cylindrical symmetry. 



idR [ d^x {xf'xfi)'+''/^ + Bd Rm f d^x {xy {x'P x^f^ + 

Joip{p,q) Jap{p,q) 

+ Bd{Y.Rkk) 



(6.55) 



/ * 

Japip,q) 



/ d'^x {xy {x'P Xpf/^ -{d-l) j d'^x {x~P Xpfl'^ . (6.56) 



This can be used to get rid of the integral involving {x ) , and get 



A2 = T 

where 
and 



2d+2 



-^((^d + ■j^^Hd,2,d/2) + ■j^(dBdHd,2,d/2 - BdHdfi,i+d/2) 



oip{j>,q) 



d'^x^x^ {xi^xay 



Hd,i,j+l/2 



= -J d''x{xy{x^x^y+<'/''- 



(6.57) 
(6.58) 
(6.59) 



here j can be either an integer or a half-integer. 

Let us now compute these coefficients. In the calculations that follows we will treat the above 
integrands simply as functions, and "forget" that their source is a curvature. If we imagine slicing 
ap{p, q) into t = constant balls and then slicing each ball into r = constant spheres, it is easy 
to see that 



/ d^x {x^ - py ((x^ - pi") {Xi, - Pi,)y/'d''x 

''oip(p,q) 



(6.60) 



r((d-l)/2) 
By changing variables to 

the above expression becomes 



( dt f f dr (f - r^l^ + f dt f f * dr r''-^(f - ry''') . 
^Jo Jo Jt/2 Jo ' 



_ i _ r 

r i 



(6.61) 



/ d<^x(xQ-p6)^((x^-/)(x^-p^)) 

Jap(p,q) 



j72 



r((d-i)/2) 



du (u^+^+'^-i dS ^-''{X - 5^)^/2) + 

+ / du(^z^+^+'^-M dss<^-2(i-s2)^/2) 



(6.62) 



Since the limits of integration in the first term are constants, that term can be represented as a 
product of two separate integrals. After evaluating the n-integral, the expression becomes 



Jap(p,q) 

27r(''-i)/2 



i/2 



m-i)/2) 



J+3+d 



_{i + j + d) 2i+o+d 



"l/u-l 



(6.63) 



+ 



J \ I'll 



From now on the calculation splits into four cases: even and odd d and even and odd j. From 
the original intentions of the calculation it is clear that whenever d is odd j is also add, and 
viceversa. So only these two cases need be considered. 

Case 1: Both d and j are even. Since j is even, denote it as 



j = 2h. 

Expanding (1 — s^)'* binomially, the integral becomes 

f d<^x {x^ - p^y ((x^ - p^) {xf, - pp)y^'' 

J a ' 



(6.64) 



oip{v,q) 

m-i)/2) 



T 



i+j+d 



fe=0 



kj [(i + 2/j + d)2'+2'»+<^ 



/' 

Jo 



dss' 



+ 



(6.65) 



+ [ duf«^+2'^+<^-i r 

Jl/2 ^ Jo 



After evaluating the s integrals this becomes 

d'^x(xO-/)^((x^-/)(x^-p^))^/' 



L 



27r('^-i)/2 



r((d-i)/2) 



r 



i+j+d 



fc=0 



+ 



+ 



2i+2h+d (i + 2/i + d) (d - 1 + 2k) 



j du (if 

Jl/2 ^ 



1/2 V d-l + 2k 

After pulling out d — 1 + 2A; and expanding out (1/n — i)d-i+2k^ ^^idX becomes 

d'^x (xO - py {{xP - p~n {x-^ - Pi,))'''' 

2-(-/-l)/2 



(6.66) 



r((d-i)/2) 



^ d-1 



k=0 



+ 2k\kJ\ 2i+'^h+d (i + 2h + d) 



+ 



(6.67) 



^ 'v^^ ,,i(d-l + 2k\ 



,+2h-2k+l 



Evaluating the integral that is left this becomes 

d-^x (x" - py ((x^ - p^) (x^ - pfi)y'^ = Haih r'^'^" , 



(6.68) 



where 



d-l+2k . ^ , _ /lNl+2h-2fe+;+l 



2k + 1 + 1 



Case 2: Both d and j are odd. Since j is odd, it will be replaced with 

j = 2/1 + 1. (6.70) 
Then the binomial expansion tells us that 

{i-sy=v^^i:(%-ifs-\ (6.71) 

fc=0 ^ ^ 

which means that the original integral can be rewritten as 



h 

X 



X 



•/ 1 /2 



(6.72) 



'1/2 

Since d is odd, so is d — 2 + 2A;. Thus, 

d - 2 + 2A; = 2a + 1 , (6.73) 

where 

^_d-3 + 2*^ (6.74) 



Thus, the integral of interest is J s^a+i vT^-7ds. 
By using 

s = sin^ (6.75) 

the integral becomes 

J Vl - ds = y sin2«+i ^ cos^d(sin^) = ^ sin^^+i cos^dS. (6.76) 
Combining one of the sin 6 factors with d^ gives 

Expanding sin^" ^ as 

sin^"^ = (1 - cos^^)" = ^i-^)^(f) (6.78) 

6=0 ^ 



the above integral becomes 



J ds = - ^ (-1)^ Q J d(cos^) cos2''+2 ^ (6.79) 



26 + 3 



By substituting s = sin ^ this becomes 

A {-If ( 



I d-s s^^^^ = - E ^ (5) (1 - ^'r'" ■ (6-81) 

Substituting 2a + 1 = d — 2 + 2A; we obtain 



I d-s ^-'^'^ = - E ^ (1 - ^"f^"" ■ (6-82) 



(-1)^ A + 

6=0 

Thus, the original integral becomes 

,?i+l/2 



^ 27r('^-i)/^ m'tT(-i)Yfc + ^\ 

r((d-l)/2)£-J^ ^ \k) 2h + A b ) 

((I __L)- _,)]}., .S3) 



{i + 2h+l + d) 2^+2/^+1+'^ 7i/2 

By evaluating the (—1) term of the integral over u, and also pulling u out of the denominators, 
the expression becomes 



f d^x (x^ - p^y ((x^ - pfi) {xi, - pf,)) 

Jap{p,q) 

t^^ \kj 26 + 3V h y V (i + 2/1 + 1 + d) 2^+2'^+i+<i 



r((d-i)/2) 



1 _ (l^«+2fe+d+l /•! 



i + 2/1 + d + 1 7i/2 
By expanding (2u — 1)''+^ we get 



- + y duir+2^+'^-2''-3(2u- 1)^+3/2)] . (6.84) 



c=0 

By setting 



dlI^Z^+2ft+d-26-3(2^_l)H3/2 

7l/2 

V(-l)'^f^ + ^')2''-^+i /' dIZlI*+2'*+'^-^-'=-2x/2¥^. 



(6.85) 



'1/2 

= V2u- 1 (6.86) 



this becomes 

du u'+^h+d-2b-^ _ 1)6+3/2 (g_87) 



1/2 



C=0 



lh+l\ r. ^^ r ^ /v^ + l\i+2h+d-b--2-c 



We can carry out the v integration if we expand the power of the binomial (u^ + 1), and we get 

du ir+2fe+<i-26-3 (2 ^ _ 1)6+3/2 89) 



c=0 ^ ^ ^ e=0 

Substituting these into the original integral we obtain 

d'^^ (x5 - p^r ((^^ - (a^M - Pm)) = ^^d,i,/.+i/2 r^'^+^+'^+i , (6.90) 

lap{j>,q) 

where 

r n.\ fc+('i-3)/2 . ..6 /r,^ <i-3^ 



/ 

Jar 



+ (6.91) 



fe=0 ^ 6=0 

1 1 _ (■l)i+2/i+<i+l 

{i + 2h + d + l) 2^+2'»+rf+i ~ i + 2/i + d + 1 



V - ^ V e / 2e + 3 

c=0 ^ ^ e=0 ^ ^ 

Finally, let's compute A3. Writing the metric in normal coordinates, 

(jfiv — Vp-i' 3 ^p,pva X 1 (6.92) 

we find that 

V|det£?| = l-lRpaxPx^ . (6.93) 

Thus, 

Asd = -jRpaJ d''xX^xmx^-p''){Xa-Pa)f\ (6.94) 

As was done with A2, the correction to the correction term will be neglected, which means that 
integration is performed over ap{p, q) instead of a{p, q) and no correction term is introduced to 

the V{a{p, x)) when the expression ((x° — p") (x,-^ — Pa))'^ was used. 
Substituting the above into the expression for A;^^^ gives 

Asd='^Rpa [ d''xxi'xmx''-p''){Xa-Pa)f\ (6.95) 

O Jap(p,q) 

By cylindrical symmetry, (x^)"^ in the above integral can be replaced with (x^)^ which means 

kd 

-"00 

'op(p>5) 
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i?00 / d^X {X^f {{X^ - p") {Xa - Pa))"''' + 



- (jlR-kk) / ((^'^ -K) (X^-P^)) 

Jap{p,q) 



d/2 



(6.96) 



Again, by cylindrical symmetry, 

Jap{p,q) Jap{p,q) 

-{d-1) [ d''x{x^f {(x^-p^) iXa-Pa)f' , (6.97) 



which allows us to express the integral involving (z^)^ in terms of integrals involving (x'^)^. 
Substituting this expression into the expression for A^d and doing some simple algebra gives 

^3d = ~ 5(j_ 1) ''"^ "^^ {Hd,o,i+d/2 - Hd,2,d/2) + Roo idHdxd/2 - Hdfi,i+d/2)\ ■ (6.98) 
Adding Ai,^, /S.2d and A'^d together, the total correction becomes 

Ad = Aid + + Am = r''^+' (Q R + Dd Roo) , (6.99) 

where 

d+2 



(6.100) 



fe=0 

and 

d+2 



(6.101) 

Now, as we were computing the corrections to the integral, we almost forgot the main term! 
Here it is: 

[ d'^x kd ((x^ - p^ {Xfi - pp)f^ = kd t'"' Hd,o,d/2 ■ (6.102) 
This means that the total integral is 

f d'^x V{a{p, x)) = kd Hd,o,d/2 r"^ + (Q i? + Dd Roo) r"^+' • (6.103) 

•Ja{p,q) 

At the beginning of this section it was shown that the gravitational pre-Lagrangian is given by 
£{^,E,p,q) = -^(v\a{p,q))+E [ d'^r ^(-1)^-^ detgV{p,r)) . (6.104) 

By substituting the above expression for the integral, as well as 

F(a(p, q)) = kdr'^ + r'^^^ {Ad R + Bd Roo) , (6.105) 

we obtain 

C{^, E,p, q) = T^'' kd {kd + E Hd,o,d/2) (6.106) 
+ T^'^+^ {R (2 kdAd + E Cd) + Roo (2 kdBd + E Dd)) . 



Thus, in order to minimize variations, we have to get rid of the 0{t ) contribution, which comes 
form the Rqq term. Thus, is selected in such a way that the coefficient of Rqo vanishes: 



2kdBd + EdDd = 0^ Ed = -2 kd Bd/Dd . (6.107) 
Finahy, substituting this into the expression for the gravitational Lagrangian gives 



AttG 



V ~d7 — JJ^ 



(6.108) 



As long as d is fixed, the first term in this Lagrangian is constant, and in fact has the form of 
a cosmological constant; its role, however, needs to be studied in more detail, since it appears 
to have a very large value, while a different type of argument based on causal sets by R Sorkin 
[9] leads to a value of the right order of magnitude for the cosmological constant. The second 
term in (I6.108P is proportional to the scalar curvature R, as expected, and can be rewritten in 
the more familiar form 

/:eh = ^^^, (6.109) 
where vq is the volume taken up by a single point of a causal set and Gd is given by 

Gd = G ^ ^^'^ . (6.110) 

2kdTl''^\AdDd-BdCd) 

Regarding the use of the above gravitational Lagrangian, it should be noted that variations 
in the partial order ^ are much more general than the variations in the metric that lead to 
the Einstein equation in the continuum version of the theory. One can speculate that, while 
causal-set variations that stay within the class of manifoldlike ones will still lead to an effective 
continuum gravitational field equation, based on the second term of (j6.108|) . the first term in 
that expression suggests that other types of variations may lead to dynamical restrictions on 
how manifoldlike "classical causal sets" are allowed to be, and particularly to conditions on their 
dimensionality. Such topics are beyond the scope of this contribution. 



7. Conclusion 

In this work, Lagrangians for the main types of bosonic fields [10] have been successfully 
defined for arbitrary causal sets. The machinery we introduced allowed us to derive otherwise 
complicated-looking causal-set Lagrangians from simple Alexandrov-set based Lagrangian 
generators, and we showed that, if the relevant fields are defined on a Lorentzian manifold and the 
Alexandrov sets are small enough that we can use the appropriate slow- variation approximation 
for the fields inside them, the Lagrangian densities coincide with the ones used in standard field 
theory. For similar work on fermionic fields, see reference [TT]. We should point out here that, 
in light of the non-linearity of the procedure for selecting an Alexandrov set, while Lagrangians 
for different fields can be added, the same is not true for Lagrangian generators. Thus, for a set 
of fields with individual Lagrangian generators Ji, we can formally write 

J = {Ju...,Jn}. Lj = Y,Lj^. (7.1) 

This line of research inherits the same problem that other approaches to quantum field 
theory have to face: if we treat gravity like all other fields, then a path integral over it implies 
an integration over all possible topologies and geometries, thus leaving us with none of the 
background information that is needed to define propagators or solve any other problems in 
ordinary quantum field theory. While this issue has been addressed in other approaches [12] . 



our belief is that in our context these should be addressed by an appropriate theory of quantum 
measurement that accommodates general relativistic covariance and includes gravity as one of 
the "measured" fields. A formulation of that theory is beyond the scope of this work. 

Another difficulty is that even if the topology was not an issue (say, we are dealing with a 
toy model of propagation of non-gravitational fields in fixed gravitational background), there 
are simply too many degrees of freedom to integrate over. In the case of regular quantum field 
theory this issue is addressed by visualizing spacetime as a regular cubic lattice, which allows 
one to perform the integrals over "all lattice points at once". In the case of a general causal 
set, however, such structure is absent, which means that one has to separately integrate over 
the fields at every single point in spacetime, one by one. While numerical simulations on the 
computer might make it possible to do so in toy models of causal sets consisting of only a few 
points, this becomes a problem for more realistic models. There is a way to make oneself feel 
better about it: The Einstein equation does not have many exact solutions either. Nevertheless, 
in order to bring causal set theory to the level of the Einstein equation, there have to be some 
simple cases for which there is a solution, such as for example a field with a very large gradient 
which assures a specific selection of Alexandrov sets. This, however, is still work for a future. 

The above problems are made much worse by the fact that causal set theory claims to be a 
candidate for a quantum theory of gravity, which implies that it should be precise enough to 
talk about the gravitational field of elementary particles: a precision much higher than the one 
of any other existing theory. Thus, reducing it to crude approximations might well defeat the 
very purpose causal set theory is there in the first place. 

Finally, the fact that the selection of Alexandrov sets is non-linear implies that the 
superposition principle does not work on microscopic scale. This might also put under question 
other principles that are based on linearity, such as the law of conservation of energy. In the 
case of smooth differentiable behavior, since these Lagrangians do reduce to the ones observed 
in regular quantum field theory, all the known laws, including the law of conservation of energy, 
should hold. They can, however, be violated once the behavior of fields is no longer smooth 
(such as, for instance, near the big bang). On the positive side, the existence of a source of 
non-conservation might have verifiable cosmological consequences. 
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